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ABSTRACT 

Strange stars are one of the possible compact stellar objects that can be formed after a supernova 
collapse. We consider a model of strange star having an inner core in the color-flavor locked phase 
surmounted by a crystalline color superconducting layer. These two phases constitute the quarksphere, 
which we assume to be the largest and heaviest part of the strange star. The next layer consists of standard 
nuclear matter forming a ionic crust, hovering on the top of the quarksphere and prevented from falling 
by a strong dipolar electric field. The dipolar electric field arises because quark matter is confined in 
the quarksphere by the strong interaction, but electrons can leak outside forming a few hundreds Fermi 
thick electron layer separating the ionic crust from the underlying quark matter. The ionic matter and 
the crystalline color superconducting matter constitute two electromagnetically coupled crust layers. We 
study the torsional oscillations of these two layers. Remarkably, we find that if a fraction larger than 10~ 4 
of the energy of a Vela-like glitch is conveyed to a torsional oscillation, the ionic crust will likely break. 
The reason is that the very rigid and heavy crystalline color superconducting crust layer will absorb only 
a small fraction of the glitch energy, leading to a large amplitude torsional oscillation of the ionic crust. 

Subject headings: stars: neutron - stars: oscillations 


1. Introduction 


The properties of hadronic matter at densities larger 
than the nuclear saturation density are mostly un¬ 
known. Even though much effort has been devoted 
to the study of very dense hadronic matter, we still 
do not know the actual ground state of matter as a 
function of the baryonic density. In particular, we are 
still unable to prove or fa lsify the hypothesis advanced 
by Bodmer and Witten ( Bodmer 1971 : Wittenl 1984 1 
that standard nuclei are not the ground state of matter. 
According to this hypothesis the energetically favored 
ground state of baryonic matter could be a collapsed 
state : an hadronic configuration corresponding to a 
(not yet experimentally observed) short-range free- 
energy minimum of the strong interaction. For small 
baryonic numbers, the collapsed state can be thought 
as a droplet of quarks and gluons having the size of 


few Fermi, smaller than the ten Fermi size character¬ 
istic of standard nuclei. The droplet has a surface ten¬ 
sion given by the bag pressure, which keeps the matter 
density inside the droplet at almost a constant value 
larger than the saturation density of standard nuclear 
matter. Therefore, with increasing baryonic number, 
A, the size of the droplet grows with r ~ A 1 / 3 , as 
characteristic of self-bound objects. In the following 
we will assume that for any baryonic density the quark 
droplet corresponds to the free-energy minimum of 
the system, meaning that it is not possible to mini¬ 
mize the free energy by fission processes. By contrast, 
standard nuclear matter can only form small clumps 
corresponding to fission stable nuclei. 

At vanishing temperature and considering that the 
only effect of the strong interaction is a bag pressure, 
the collapsed state can be thought as consisting of al- 
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most free quarks filling the pertinent energy levels up 
to the corresponding Fermi energy. When the light 
quark chemical potential exceeds the strange quark 
mass, the strange quark states start to be populated by 
means of weak decay processes. In this case, the col¬ 
lapsed state co rresponds to “catalyzed” u, d, s quark 
matter ( Witten 198 4). the so-called strange matter. A 
very big clump of strange matter is called a strange 
star ( Alcock et al.lll98(l Haensel et aL 1986 1. Strange 
stars and neutron stars can be viewed as two different 
end products of supernova explosions and are classi¬ 
fied as compact stellar objects (CSOs), which are ob¬ 
served as stars having the radius of about 10 km and 
the mass of about a solar mass, Mq. Actually, assum¬ 
ing that strange matter is self-bound implies that the 
size of a strange star has no lower bound, meaning that 
strange stars much smaller than standard neutron stars 
can exist. Gravity plays a role only for very massive 
objects, restricting the mass to about 2 Mq (for a suf- 
fici ently stiff equation of s tate (EoS)), see for exam¬ 
ple Mannarelli et al.l ( 20141) . 

Assuming that the ground state of hadronic mat¬ 
ter consists of strange matter, it still does not clar¬ 
ify unambiguously the properties of the system. The 
largest value of the quark chemical potential that can 
be reached in massive strange stars is of the order 
of 400 — 500 MeV. Even considering this extreme 
case, the strong interaction is still nonperturbative and 
quantum chromodynamics (QCD) is not under quan¬ 
titative control. Therefore, approximation schemes 
must be used. Analyses using various models indi¬ 
cate that at the densities relevant for strange stars, de- 
confined and cold quark matter isjikely in acolor su¬ 
perconducting (CSC) phase , see Raiagonal & Wilczekl 
( 2000l> : Alford et al. ( 2008 ): Anglani et al. ( 2014 ) for 
reviews, in which quarks form Cooper pairs breaking 
the SU (3) C oior gauge symmetry. The reason is that the 
estimated critical temperature of color superconduc¬ 
tors is at least of the order of few MeV, much larger 
than the tens of KeV temperature of few seconds old 
CSOs. Thus, if quark matter is present in CSOs, it 
should be in a CSC phase. 

The CSC phase is actually a collection of phases 
and pinning down the favored quark pairing is not 
trivial. Using models based on one gluon exchanges 
or on instanton exchanges, it can be shown that at 
asympt otic densities the color-flavor locked (CFL) 
phase dAlford et al.l 119991) is energetically favored. 
In the CFL phase u , d, s quarks of all colors pair 
coherently maximizing the free-energy gain. Al¬ 


though this phase is very robust, it might be that 
strange stars do not completely consist of CFL mat¬ 
ter. The reason is that when the effective strange quark 
mass has a value comparable with the quark chemi¬ 
cal potential, a considerable free-energy penalty re¬ 
sults for producing strange quarks. If the free-energy 
penalty is larger than the free-energy gain associated 
to the CFL pairing —exceedi ng the corresponding 
Chandrasekhar- C logston limit (IChandrasekhan 1 1962 ; 


Clogston 1962 ; Anglani et al. 20141) — a different 


phase is favored. In the present paper we assume that 
the CFL phase is realized in the central and denser 
part of the strange star and the next favored phase 
down in density is the crystalline color superconduct 


ing (C CSC) phase ( Alford et al. 2001 at Anglani et al 


201 4), which is then realized in the outer and less 


dense part of the strange s tar dRupak & Jaikumar 


201 2: Mannarelli et al.l 2014 1). The CCSC phase is 


characterized by a periodic modulation of the diquark 
pairing. For our purposes, the most important property 
of the CCSC phase is that this periodic modulation 
is mecha nically rigid, with an extremely large shea r 
modulus ( Mannarelli et al.il2007t Anglani et al. 2014 ). 
In our model the CFL phase and the CCSC phase form 
the quarkspliere comprising most of the mass of the 
strange star. It is known by various mo del calculations, 
see for example Anglani et al. ( 20141) . that the CCSC 
phase is the favored phase for sufficiently mismatched 
quark Fermi spheres, corresponding to a quark chem¬ 
ical potential in a certain range of values. The actual 
matter density at which the CCSC is favored depends 
on the detailed dependence of the strange quark mass 
on the quark chemical potential and the corresponding 
value of the CFL gap parameter. Since these quantities 
cannot be precisely computed, we will assume that at a 
certain radius, Rqfl, there exists a boundary between 
the CFL and the CCSC phase. 

An important aspect of dense quark matter in CSOs 
is that it must be electrically neutral. The typical com¬ 
bined effect of charge neutrality, j3 decays and non¬ 
vanishing strange quark mass, M s , is t o populate the 
electron states, see for example Alford et al.l (12001b ). 
Qualitatively, the reason is that a large strange quark 
mass disfavors the appearance of strange quarks by 
light quark j3 decay, thus the electrical neutrality con¬ 
dition is satisfied by populating electron states. A no¬ 
table exception is the CFL phase, in which the sym¬ 
metric pairing induced by the strong interaction forces 
equal number of u, d and s quarks, thus no electrons 
are present. On the other hand, in the CCSC phase, the 
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quark Fermi momenta are mismatched and electrons 
are needed to maintain the electric charge neutrality. 
Therefore, our model of strange star consists of a cen¬ 
tral electron-free region of CFL matter in contact with 
an electron-rich region of CCSC matter. The presence 
of electrons in the CCSC phase leads to an interesting 
phenomenon happening at the surface of the strange 
star where the strong interaction confines quarks be¬ 
neath the CCSC strange star surface. Electrons are 
only bound by the electromagnetic force and will typi¬ 
cally spread outside the quarksphere for a length-scale 
of the order of the Debye screening length, of hun¬ 
dreds of Fermi. A CSO made by quark matter sur¬ 
rounded by an electrosphere is called a bare strange 


Alcock et al. ( 19861) for more details. 


star, see 

In some circumstances it is possible that a standard 
ionic crust hovers on the top of the strange star. The 
reason is that in the formation process, or by accretion, 
the strange star attracts hadronic matter. Neutrons are 
absorbed by the quarksphere, but ions are repelled by 
the positively charged quark surface. Thus, ions hover 
on the top of the strange star as far as the repulsive 
electric force and the gravitational force balance. If 
the accreted material is enough, it will form a ionic 
crust separated by the deconfined quark matter surface 
by a thin electron layer. This CSO is called a nonbare 
strange star (lAlcock et al.lll986l) . In the following we 
will consider a model of nonbare strange star having 
two crusts, an inner crust layer consisting of CCSC 
matter and an outer crust layer consisting of standard 
ionic matter. 

In lMannarelli et al.l (120141) we studied the torsional 
oscillations of bare strange stars. In the present pa¬ 
per we extend that study to nonbare strange stars, an¬ 
alyzing the torsional oscillations of the two coupled 
crust layers. An intriguing aspect of this model is that 
these two crusts have very different shear moduli and 
densities, therefore nontrivial dynamical process can 
happen. We focus our analysis on £ = 1 modes, cor¬ 
responding to oscillatory twists of the crust. These 
modes do not conserve angular momentum, therefore 
we assume that they are triggered by events that trans¬ 
fer angular momentum to the strange star crust. A typi¬ 
cal event of this sort is a stellar glitch. For definiteness, 
we will assume that an energy of the order of the one 
released in a Vela-like glitch is conveyed to the t =1 
modes; we will show how our results can be appropri¬ 
ately rescaled if a different energy scale is used. One 
of the most interesting results that we obtain is that the 
shear strain has a radial dependence with a maximum 


much closer to the star surface than in standard neutron 
stars. Moreover, in standard neutron stars the energy 
of the torsional oscillations is spread across the entire 
crust, which is more than a km thick. In the considered 
model of nonbare strange stars, for a sufficiently thin 
CCSC crust, all the energy of the torsional oscillation 
is conveyed in a layer few hundred meters thick, corre¬ 
sponding to the ionic crust at densities below neutron 
drip. Therefore, the ionic crust layer will likely crack 
even if a small fraction, of order 10 -4 , of the Vela-like 
glitch is conveyed to the ionic crust. 

The present paper is organized as follows. In Sec. [2] 
we present a general description of the nonbare strange 
star model. In Sec.[3]we determine the stellar structure 
by solving the hydrostatic equilibrium equations. In 
Sec. [4] we discuss the torsional oscillations of the stel¬ 
lar crust. In Sec. [5] we draw our conclusions. In Ap¬ 
pendix |A] we discuss a simple toy model of two rigid 
slabs, determining the corresponding torsional eigen- 
frequencies. 


2. General description 


The considered nonbare strange star model is de¬ 
picted in Fig. |T] The u. d and s quark matter is radially 
confined by the strong interaction in the quarksphere, 
with radius R q . This part of the star can be considered 
as a big hadron with an extremely large baryonic num¬ 
ber. Since the quark chemical potential is large and 
the temperature is sufficiently low, deconfined quark 
matter is assumed to be in a CSC phase. The inner 
part of the quarksphere, up to the radius Rqfl = aR q , 
with a < 1, is in the CFL phase. The CFL phase is sur¬ 
mounted by the CCSC phase, extending between Rqfl 
and R q . The outer part of the nonbare strange star con¬ 
sists of a thin electron layer, extending between R q 
and R P . This region is on ly a few hundreds of Fermi 
thick (lAlcock et al.l 19861) , thus when considering the 
hydrostatic equilibrium configuration we will not dis¬ 
tinguish between R q and R e . 

The next layer on the top of the electrosphere is the 
ionic crust, which is composed by ions and electrons. 
This layer extends betwee n R P and R, reachin g at most 
the neutron drip density ( Alcock et al. 19861) . Densi¬ 


ties above the neutron drip point cannot be attained 
because the ionic crust is maintained from collapsing 
on the underlying strange star surface by the electro¬ 
static surface field. If the density of the crust reaches 
the neutron drip point, neutrons are liberated falling 
on the underlying strange star surface by gravitational 
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Fig. 1.— Pictorial description of the strange star struc¬ 
ture (not in scale). The stellar radius of the considered 
model is R ~ 9.18 km. Most of the strange star, for 
0 < R < R q ~ 8.98 km, consists of u,d,s quark 
matter in a color superconducting phase, forming the 
quarksphere. The inner part of the quarksphere, for 
0 < R < Rc fl, consists of electron-free color-flavor 
locked u,d,s matter; the outer part of the quarksphere, 
Rcfl < R < R q , consists of the electron-rich crys¬ 
talline color superconducting u, d. s matter. The radius 
Rcfl is unknown and we consider values in the range 
0 < -Rcfl < Rq- On the top of the quarksphere the 
electron layer extends between R q and R e ; it is just 
few hundreds of Fermi thick and consists of the elec¬ 
trons leaked outside the quarksphere. The last layer is 
the ionic crust, which has the same properties of stan¬ 
dard neutron star crust. 


attraction and are thus absorbed in the quarksphere. 
This process pu ts a l imit o n the maximum mass of 


the ionic crust ( Alcock et al. 19861) . The ionic crust 


layer of nonbare strange stars has the same mechan¬ 
ical properties of the outer crust of standard neutron 
stars, meaning that we will considered the known EoS 
of nuclear matter below neutron drip and the estimated 
values of the nuclear matter elastic moduli. Moreover, 
as we will discuss below, the ionic crust is in part liq¬ 
uid forming the so-called ocean. 


One important aspect to clarify is that for the back¬ 
ground configuration we assume hydrostatic equilib¬ 
rium, but oscillations are described considering an 


elastic response to the mechanical stress. Although 
the descriptions of the equilibrium and of the mechan¬ 
ically excited state of the star seem in contradiction, 
they can be justified as follows. The hydrostatic equi¬ 
librium is appropriate if a liquid description holds. 
This is a reasonable approximation in any layer of 
the strange star. As is well known, for matter den¬ 
sity above 10 5 g/cm 3 the pressure in the ionic crust is 
mainly due to degenerate electrons that can be treated 
as a fluid. The remaining part of the ionic crust, hav¬ 
ing smaller density, has a very small mass that can be 
safely neglected in the solution of the hydrostatic equi¬ 
librium equations. Regarding quark matter, the CFL 
phase is a fluid and thus a liquid description is appro¬ 
priate. The CCSC phase is rigid because of the rigidity 
of the periodic pattern of the gap parameter, however 
the equilibrium pressure of quark matter does weakly 
depend on the presence of the condensate. We shall 
assume that such a contribution can be absorbed in a 
liquid-like description of quark matter. 

Regarding the mechanical response, we will only 
include the effect of the shear modulus, because we 
focus on torsional oscillations that only have a tangen¬ 
tial stress component. Obviously, these oscillations are 
confined to the crusts of the nonbare strange star, be¬ 
cause a fluid cannot support shear waves. The vari¬ 
ous interfaces between the different layers will be de¬ 
scribed by appropriate boundary conditions. Assum¬ 
ing that a small shear perturbation acts on the system, 
we expand the displacement field, the pressure and the 
matter density respectively as follows 

U = U 0 +u P ik = P 0 5 ik -Iii k p = p Q +Sp,(l) 

where the quantities with a subscript 0 correspond to 
the time independent background and the remaining 
quantities represent the linear perturbations. 


3. Background configuration 


We set as the equilibrium configuration the one with 
Uq = 0 and assume negligible background magnetic 
field. Therefore, our analysis is strictly valid for slowly 
rotating nonmagnetar stars; however we will briefly 
comment on the effect of a nonvanishing magnetic 
field. In the stationary state, electrons and quarks are 
confined in regions with a net charge and are in ther¬ 
mal equilibrium. This is possible because nonneutral 
systems in appropriate geometries can be confined and 
in thermodynamic equili brium; a typical example are 
non-neutral plasmas, see Dubin & O’Neil ( 19991 ) for 
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Fig. 2.— Density profile of the considered model of 
nonbare strange star. Strange matter is present for 
r < R q ~ 8.98 km; the solid (blue) line represents 
the quark matter density. The ionic crust layer extends 
for R q < r < R, with R = 9.18 km the star radius; 
the dashed (black) line represents the nuclear matter 
density. 


a review, which, unlike quasineutral plasmas, can be 
brought to equilibrium and confined. 

The equilibrium is determined by the balance be¬ 
tween the hydrostatic pressure and the gravitational 
attraction, which is appropriately described by the 
Tolman-Oppenheimer-Volkov (TOV) equation 


dp G(p + p) (to + 47t p r 3 ) 
dr r(r — 2 G m) 


( 2 ) 


where m(r) = dr' r' 2 p(r') and G is the gravi¬ 
tational constant. Two different EoSs must be used, 
Pqm{p), for quark matter in the quarksphere and 
Pnm(p), for nuclear matter in the ionic crust. In both 
cases we assume that the temperature is so low that 
the background distribution can be approximated by a 
Fermi liquid at zero temperature. Quarks are strongly 
interacting and a first principle calculation of the EoS 
is unfeasible. However, if the leading contribution of 
the quark interaction is quark pairing, we expect cor¬ 
rections of the order A 2 / p 2 < 10% to the free Fermi 
gas EoS. To effectively take into account the strong 
interaction we use the general parame ter izati on of the 
quark matter EoS given in Alford et ak (2005) 


Dqm — 


47T 2 


<24/4 -f- + -®eff ! (3) 


where 04 , a 2 and B e g are independent of the average 
quark chemical potential //. We use the set of param¬ 
eters 04 = 0.7, 02 = (200 MeV ) 2 and f? e ff = (165 
MeV ) 4 (similar r esults hold for the two se ts of param¬ 
eters discussed in lMannarelli et al.1 (120141) ). 


For the ionic crust EoS we assume that it con¬ 
sists of a Coulomb crystal embedded in a degen¬ 
erate electron gas. We only need an expression 



the highest density of nuclear matter corresponds to 
the neutron drip point. This assumption corresponds 
to the maximization of the ionic crust radius. More 
refined studies show that properly taking into account 
the electronic pressure results in a reduction of the 
ionic crust mass and thickness (lMartemvanovT ll994: 


Huang & Lu 1997). However, the detailed extension 


of the ionic crust is not relevant for our purposes, we 
are only interested in order of magnitude estimates of 
the effect of a ionic crust. Therefore, in our approach 
the radius, R q , corresponding to the boundary between 
the quarksphere and the ionic crust, is determined by 
the condition that the pressure of quark matter equals 
the pressure of nuclear matter at the neutron drip point, 
i.e. PQM(Rq) = PNM(Rq) = PND — 7.8 X 10 29 
dynes/cm 2 . The star radius, R, is determined by the 
boundary condition on the pressure pnm(R) = 0 , 
meaning that the pressure at the surface of the ionic 
crust vanishes. Note that the effect of the magnetic 
field on the matter distribution is small even consid¬ 
ering the extreme values characte ristics of mag n etars, 
see for example the discussion in Frieben & Rezzollal 

d2012h. 


In Fig. U we report the matter density profile for 
the considered strange star model obtained solving the 
TOV equation for a strange star having mass of about 
1.4 Mg. The obtained star radius is R ~ 9.18 km, 
corresponding to the outer surface of the ionic crust 
layer. The matter density in the quarksphere is roughly 
constant, changing by less than a factor 2. There is a 
jump at i? = ~ 8.98 km corresponding to the tran¬ 

sition between quark matter and the ionic crust. The 
nuclear density changes by about 10 orders of mag¬ 
nitude across the ionic crust layer, from the iron-like 
density of ppe — 7.8 g/cm 3 to the neutron drip density 
Pnd — 4.3 x 10 11 g/cm 3 . This is the same radial de¬ 
pendence found in the outer crust of standard neutron 
stars. 
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4. Torsional Oscillations 

We now discuss the fluctuations of the displacement 
fields. Assuming a fluid description, the oscillations 
obey the Euler’s equation 

P -^ + d k n ik = 0 , (4) 

at 

and the continuity equation 

d t Sp+ V- (pu) = 0. (5) 

For simplicity we have not considered dissipative pro¬ 
cesses and assumed that the matter density of the vari¬ 
ous species does not change, in other words the system 
is in “chemical” equilibrium. We have also omitted 
the collision term, meaning that we assume that colli¬ 
sions are so fast that local equilibrium is reached in a 
time scale negligible compared with the one of exter¬ 
nal forces. All the perturbations oscillate at the same 
frequency <r; in particular we focus on the eigenmodes 

u = e i ” t £, (6) 


where Vi is the shear velocity and u, is the shear mod¬ 
ulus. The index i = 1,2 refer to the CCSC crust layer 
and to the ionic crust layer, respectively. The differ¬ 
ential equations describing the oscillations of the two 
crusts are decoupled, because no long range forces are 
present. However, a coupling between the two oscil¬ 
lations is produced by the boundary conditions (BCs), 
effectively describing the effect of short-range electro¬ 
magnetic forces between the two crusts. The set of 
BCs that we use for solving the two Euler’s equations 
is the following 







W 1 (R q ) = W 2 (R q ) 



( 10 ) 



-Rq 

( 11 ) 

( 12 ) 


(13) 


and assume that all of the fluctuations are proportional 
to the displacement field, meaning that any perturba¬ 
tion vanishes for £ = 0. We restrict the analysis to 
torsional oscillations, corresponding to transverse os¬ 
cillations with no radial displacement 


V • u = 0 


U r = 0. 


(7) 


Upon plugging the above expressions in Eq. (0 it can 
be easily shown that no matter density fluctuation is 
produced. 

For nonrotating and nonmagnetic stellar models the 
torsional oscillations do not couple with any other star 
oscillation, thus an eigenmode analysis is possible. 
According with McDermott et alJ (1 1 9881) we will in¬ 
dicate with pt n the torsional mode having harmonic 
index t and n nodes. The f-mode displacement field 
can be decomposed as follows 


Ze = W 


1 dY lm 
sin 0 d(j) 


U = -w 


dYi m 

d6 


( 8 ) 


K ?! 

(H983f) : IAndersson et al.l (12002ft for general relativistic 
discussions) the Euler’s equation in spherical coordi¬ 
nates reads 


where f? 2 — 9.15 km corresponds to the bound¬ 
ary between the solid crust and the ocean (at p ~ 
10'g/cm 3 ). These equations describe the no-traction 
BCs, Eqs. dm CD , ©. and the no-slip BC, 
Eq. dE}- The no-traction BCs simply mean that no 
net force acts on any boundary between two adjacent 
layers. The no-slip BC means that the displacement 
of the two crust layers at the interface is the same. 
Thus, we are assuming that the static friction at this 
interface is so strong that at any time the two crust 
oscillations have the same amplitude and frequency, 
a. Note that the assumption that the oscillations in 
the two crusts have the same phase and frequency was 
implicitly done in Eq. ©. 

At the CCSC-CFL boundary we only impose the 
no-traction BC, Eq. (flOt . but the two materials have 
free slip. The reason is that the CFL phase is a neutral 
superfluid, therefore it should not react with a force to 
the oscillation of the CCSC internal surface. In this 
case we are assuming that no vortices are created by 


and in the Newt onian ap proxim ation (see Schumaker & Thorne the oscillation of the CCSC boundary and that the tem¬ 
perature is so low that thermal fluctuations in the CFL 
superfluid can be neglected. A similar reasoning ap¬ 
plies to the ionic crust-ocean boundary, which is de¬ 
scribed by the no-traction BC, Eq. (IT2l) and free slip. 

The shear moduli and densities of the two crust lay¬ 
ers are fundamental ingredients for determining the so¬ 
lutions of Eq. ©. We already determined the matter 


2 W, =i 


dlogi/j, f dWi Wi 


dr 


LA. 

r 2 dr 


\ dr 
,dW, 


dr 


£(£+ 1) 


Wi 


(9) 
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density profile in Sec.0 Regarding the shear modulus 
of the CCSC phase, it can be obtained from the low- 
energy Lagrangian description in terms of phono n¬ 
like excitations, s ee Casalbuoni et alJ ( 2001 . 2002blla ): 
Mannarelli et al.l ( 20071) . In particular, the coefficients 
of the quadratic terms of the low-energy Lagrangian 
determine the linear response of the crystal to excita¬ 
tions. Unfortunately, to extract the shear modulus var¬ 
ious ap proximations have to b e used: the procedure 
used in iMannarelli et al.l (120071) relies on a Ginzburg- 
Landau (GL) expansion of the Nambu-Jona Lasinio 
model used to mimic the properties of quark mat¬ 
ter at the relevant baryonic densities. It is known 
tha t this procedure is not under quantitative control. 


see 


Mannarelli et al.l ((2006), moreover the recent anal¬ 


yses of two flavor systems dCao et al . 2015|) seem to 
indicate that many terms in the GL expansion must be 
included to have a controlled approximation scheme. 
U sing the GL expansion u p to order A 2 it was found 


Mannarelli et al. (2007) that 


z'ccsc — v o 


a (m) 

lOMeV 


) 2 (iookv) 2 ' (14) 


where 

MeV 

vo = 2.47——g- , (15) 

fm‘ 

will be our reference value. We will assume that the 
shear modulus is constant within the CCSC crust and 
hereafter we will identify the CCSC shear modulus 
with the reference value. It is certainly true that uqcsc 
depends on the quark chemical potential /i. However, 
in the interior of compact stars the quark chemical po¬ 
tential is almost constant. Moreover, as we will see 
below, the most relevant aspect is that v ccsc is much 
larger than the shear modulus of the ionic crust. Re¬ 
garding the quark matter density, it does not strongly 
depend on the radial coordinate. Therefore, we will 
consider the constant value pqm = 10 15 g/cm 3 ; a typ¬ 
ical quark matter value, see Fig. [2] 


The shear modulus of the ionic crust depends on the 
particular crystalline structure considered and on the 
plane of application of the shear stress. Early calcula¬ 
tion of the shear_modulus_of monovalent crystals were 
performed bv (Fuchs ( 1936 ). In compact stars the ori¬ 
entation of the crystals is unknown, however it is possi¬ 
ble to define an effective she ar modulus, i/ R ff, averagin g 
over directions, as shown in lStrohmaver et al.l dl99ll) . 
which should give an excellent approximate result if 
the crust has a polycrystalline structure. From dimen¬ 
sional analysis and considering that the rigidity is due 


to the electromagnetic interaction between ions, it is 
clear that zz e ff(r) oc (Z(r)e) 2 njv(r) 4 ^ 3 , where Z{r) is 
the radial dependent proton number and njv(r) is the 
number density of nuclei. The order of magnitude esti¬ 
mate of the shear modulus gives a very large value, of 
about IQ 30 dyne cm~ 2 for the inner crust, as already 
noted in Smolucho wskil dl97C)l). Using Monte Carlo 
simulations, see Strohmaver e£ak|j[l991 ), or Molecu - 
lar Dynamics methods, see iHoffman & Hevll (2012 k 
eventu ally including quantum fluctuations as in lBaiko 
(2011[), one can estimate the proportionality factor, 
finding that 


v e s(r) = c 


n N (r)(Z(r)e) 2 
a(r ) 


(16) 


where a(r) = (3/(47mAr(r))) 1 / 3 is the average inter¬ 
ion spacing and c ~ 0.1 gives an approximate result 
at any density larger than 10 6 g/cm 3 , where electrons 
form a degenerate Fermi gas and the Coulomb crystal 
model can be applied. At smaller density, the Coulomb 
crystal model cannot be applied and Eq. (fl6t does not 
apply. As an example, extrapolating Eq. (fl6l ) to the star 
surface does not reproduce the known shear modulus 
of iron. 


Summarizing, the matter densities and shear moduli 
of the nonbare strange star are approximately given by 


r 0 

for r < R c fl 


II 

0 

for 7?cfl < r < R q 

(17) 

{ VeB (r) 

for i? q < r < i? 2 


and 



f PQM 
\ Pnm(V) 

for Rqfl < r < R q 
for R q < r < R 2 

(18) 

with pq M = 10 15 g/cm 

3 and the CCSC shear modulus 


given in Eq. (fl5l ). For a realistic description of am¬ 
plitude and frequency of the torsional oscillations, the 
radial dependence of the nuclear matter density and 
of the effective shear modulus must be appropriately 
taken into account in Eq. ©• However, to disentan¬ 
gle the various aspects of the problem in Sec. 14. II we 
will approximate the ionic crust as an homogeneous 
system. We turn to a discussion of the inhomogenous 
ionic crust layer in Sec. 14.21 


4.1. System of two homogeneous crusts 

Let us approximate the nonbare strange star crust 
as made by two homogenous crust layers. We con¬ 
sider two different parameter sets for characterizing 
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the properties of the ionic crust. Parameter set A cor¬ 
responds to the nuclear matter density and the effec¬ 
tive shear modulus at the inner surface of the ionic 
crust, respectively given by Pnm(-R 9 ) = Pnd and 
v e ff(Rq) — 3.4 • 10 -5 MeV/fm 3 . The parameter set 
B corresponds to the nuclear matter density and the ef¬ 
fective shear modulus near the ocean surface, respec¬ 
tively given by p(R 2 ) = 2.8- 10 9 g/cm 3 and v e «(R 2 ) = 
2.1 • 10” 10 MeV/fm 3 . 

Since both crust layers are homogeneous, Eq. Q 
can be cast in the form of two Bessel’s equations, and 
we can determine an analytic expression for the oscil¬ 
lations for any l. In the following we focus on £ = 1, 
but similar results hold for different values of l. Im¬ 
posing the BCs given in Eq. (ITOl) and (Il3l) we find that 



Wi{r) = -^=Ci (J 3 / 2 (ra/Vi) - K z Y 3/2 (r<j/vi)) , 

(19) 

where J n and Y„ are Bessel functions of the first and 
second kind, respectively, and 

Ki = J 5 / 2 {aRia/v 1 )/Y 5/2 (aR 1 a/vi) (20) 

and 

K 2 = J 5/2 (Ra/v 2 )/Y 5/2 (Ra/v 2 ). (21) 

By the two remaining BCs, Eq. (ITO and Eq. we 
can determine the ratio C\/C 2 and the eigenmode fre¬ 
quency. 



a 


In Fig.[3]we report the dependence of the frequency 
of the first three eigenmodes on a. The solid blue line 
corresponds to the \t\ mode; the dashed red line corre¬ 
sponds to the \ t 2 mode and the dotted green line corre¬ 
sponds to the \t 3 mode. The results reported in the up¬ 
per panel are obtained with the parameter set A, while 
the results reported in the bottom panel are obtained 
with the parameter set B. 


Let us focus on the fundamental i 1 \ mode obtained 
with the parameter set A. For a < 0.75 the frequency 
of this mode increases with increasing a. Then, for 
a > 0.75 the frequency becomes almost independent 
of a. For a < 0.75 the frequency and the amplitude of 
th e oscillation are very sim ilar to the results obtained 


in IMannarelli et all (2014) for a bare strange star. In 
this range of values of a the system behaves as if the 
ionic crust layer is absent. For a > 0.75 the opposite 
behavior happens, with torsional oscillation segregated 
in the ionic crust layer with a frequency determined by 
the mechanical properties of the ionic lattice. Indeed, 
this frequency corresponds to ui 2 ~ n/2 ■ v 2 /(R — 
R q ) as described in the Appendix lAl using a simpler 


Fig. 3.— Frequency of the /-modes with t = 1 as a 
function of a for different number of nodes: n = 1 
mode solid blue, n = 2 mode dashed red, n = 3 mode 
dotted green. The plots have been obtained consid¬ 
ering both the CCSC crust layer and the ionic crust 
layer as homogeneous. For the ionic crust we consid¬ 
ered two sets of parameters. Top: v 2 ~ 3.55 x 10 8 
cm/s and p = p ND . The frequency of the 1 1 1 mode 
increases with a until a ~ 0.75. For a > 0.75 
this mode is an almost pure ionic-crust oscillation and 
does not depend on the extension of the CCSC crust 
layer. A similar phenomenon happens for the 3 t 2 
mode, which becomes an almost pure ionic-crust os¬ 
cillation for 0.5 < a < 0.75, and for the 1/3 mode, 
which becomes an almost pure ionic-crust oscillation 
for 0.85 < a < 0.9. Right panel: v 2 ~ 1.1 x 10 7 cm/s 
and p ~ 1.5 g/cm 3 . The frequencies of the modes ifi, 
1/2 and 1/3 are all very weakly dependent on a. 
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Fig. 4.— Amplitude of the oscillation of the it\ mode 
for the parameter set A for two different values of a. 
Top: results obtained for a = 0.4. Both crust layers os¬ 
cillate with a comparable amplitude. Bottom: results 
obtained for a = 0.8. The oscillation is basically seg¬ 
regated to the ionic crust. In both cases the amplitudes 
are obtained assuming that the energy of a Vela-like 
glitch is conveyed to the Ri mode oscillation. 

model with planar geometry. For a* = [1 — 2vi/R q ■ 
( R — R q )/v 2] ~ 0.75 there is a crossing between two 
different types of mode oscillations, meaning that the 
frequency of the first CCSC crust oscillation and of the 
first ionic crust oscillation are equal. Therefore, at this 
point there is a crossing between the frequency of the 
n = l mode and of the n = 2 mode. The transition 
between the two behaviors corresponds to a node at the 
interface between the two crust layers. 

The modes with more nodes have a similar behav¬ 


ior. For the it 2 mode there are two crossing points, 
one with the mode with n = 3 and one with the mode 
with n = 1. As can be seen in the upper panel of Fig. 0 
there exist two regions in which the frequency of R 2 
mode is almost independent of a and two regions in 
which it is strongly dependent on a. The it 3 mode has 
a similar behavior, as well. We have checked that for 
any considered mode, the range of values for which the 
frequency is almost constant corresponds to torsional 
oscillation basically confined to the ionic crust layer. 
Note that for any value of a there is at least one of the 
three modes that is weakly dependent on a, meaning 
that there is at least one mode that is segregated to the 
ionic crust. For the set of parameter B, the frequencies 
of the three considered modes are independent of a for 
any value of a, see the bottom panel of Fig. [3] There¬ 
fore, for the parameter set B the low lying f-mode fre¬ 
quencies are completely determined by the mechanical 
properties of the ionic crust layer. 

To better understand the difference between the 
two regions of frequencies, respectively dependent and 
constant in a, let us discuss in detail the amplitude of 
the it 1 mode. In Fig. 0]we show the amplitude of 
the it 1 mode obtained for the parameter set A, con¬ 
sidering a = 0.4 upper panel, and a = 0.8, bottom 
panel. To determine the absolute values of this ampli¬ 
tude we assume that all the energy of a Vela-like glitch, 
-EVeia ~ 5 ■ 10 42 ergs, is conveyed to the 1 1 j mode. If 
a fraction a of the Vela-like glitch is considered, then 
the amplitude has to be scaled by \fa. For a = 0.4, 
we obtain that the CCSC crust and the ionic crust os¬ 
cillate with similar amplitude that is order of tens of 
centimeters. For a = 0.8 the oscillation of the CCSC 
crust is negligible and the mode is segregated to the 
ionic crust layer. The energy of the glitch conveyed in 
the thin ionic layer generates amplitude order of tens 
of meters. This result confirms the expectation that 
t-modes having an a independent frequency are basi¬ 
cally ionic crust oscillation. We can therefore classify 
the nonbare strange star oscillations as CCSC crust os¬ 
cillations and ionic crust oscillations. For the set of 
parameters B the dependence from a of the first three 
modes is negligible and the corresponding amplitude is 
always confined in the ionic crust layer with a behavior 
very similar to the one reported in the bottom panel of 
Fig. HI However, the absolute value of the amplitude 
in this case is extremely large and can reach values 
of kilometers at the star surface. Since these values 
of the oscillation amplitude are well beyond the linear 
approximation region, let us understand the physical 
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Fig. 5.— Absolute values of the shear strain, defined in 
Eq. C3, across the ionic crust layer. Top: results ob¬ 
tained with parameter set A and for a = 0.8. Bottom: 
results obtained with parameter set B and for a = 0.4. 
In both cases the strains are obtained assuming that the 
energy of a Vela-like glitch is conveyed into a it± mode 
oscillation. 

behavior of the system considering the deformation of 
the crust. 

A measure of a solid deformation is the shear strain 
(also called the shear traction). The shear strain de¬ 
pends on the angular coordinates, but for simplicity we 
define 


which only depends on the radial coordinate; the shear 
strain at any angle can be obtained multiplying |s| by 
the appropriate angular function, see Eq. ©. We focus 


on the ifi mode and in Fig.[5]we report the correspond¬ 
ing shear strain across the ionic crust (the deformation 
of the CCSC crust layer is always much smaller). The 
upper panel of Fig. 0 corresponds to the ionic crust 
strain produced by a \t\ ionic crust oscillation with 
the parameters set A and a = 0.8, namely the re¬ 
ported strain corresponds to the amplitude shown in 
upper panel of Fig. |4] The bottom panel of Fig. |5]cor- 
responds to a ionic crust strain induced by a oscil¬ 
lation for the parameters set B, showing the extreme 
values for the deformation discussed above. Note that 
in both panels of Fig. 0 the strain has a maximum at 
the boundary between the CCSC crust layer and the 
ionic-crust layer. The shear strain is a monotonic de¬ 
creasing function of the radial coordinate because the 
ionic crust layer is homogenous and because the shear 
strain at the surface has to vanish because of the no¬ 
traction BC. As we will see in the next section, con¬ 
sidering an inhomogenous matter density results in a 
nonmonotonic shear strain. For all considered cases 
the deformation of the ionic crust layer is large mean¬ 
ing that the linear approximation is not valid. Indeed, 
in this case the crust very likely cracks before reach¬ 
ing such a large deformation. We will discuss crust 
cracking in Sec. [5] we now turn to a more realistic de¬ 
scription of the ionic crust layer. 

4.2. Inhomogeneous ionic crust 

We now include the radial dependence of the me¬ 
chanical properties of the ionic crust matter, while 
keeping both the matter density of quark matter den¬ 
sity and the shear modulus of the CCSC crust layer 
constant. In particular, we numerically solve Eq. © 
including the radial dependence of the ionic crust shear 
velocity, V 2 (r), and shear modulus, ^(r). The fre¬ 
quencies of the first three torsional eigenmodes are re¬ 
ported in Fig. [6] This figure is rather similar to the 
one reported in the upper panel of Fig. [3] showing that 
the parameter set A is the one that gives a reasonably 
good approximation of the ionic crust, as far as f-mode 
oscillation are concerned. 

Let us focus on the \t\ mode. For small values of 
a the CCSC crust oscillation frequency strongly de¬ 
pends on a. The corresponding oscillation amplitude 
is mostly confined in the CCSC crust, see the upper 
panel of Fig. [7] for a representative behavior obtained 
with a = 0.4. This amplitude is very similar to the 
one reported in the upper panel of Fig. |4j the reason 
is that the same values of the CCSC crust density and 
shear modulus have been used. Increasing the value 
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of a one reaches a critical value around a* ~ 0.7 for 
which there is a crossing between the n = 1 and the 
n = 2 modes. For a > a* the amplitude of the oscil¬ 
lations is mainly confined in the ionic crust layer and 
the oscillation frequency is almost independent of a. 
In the bottom panel of Fig. [7]is reported the amplitude 
of the iti ionic-crust oscillation obtained for a = 0.8. 
In this case, the qualitative behavior of the oscillation 
amplitude is very similar to the one reported in the bot¬ 
tom panel of Fig.[4j but the amplitude of the oscillation 
is now larger by about an order of magnitude. The rea¬ 
son is that in this case the light matter at the top of the 
ionic crust is easily displaced by the torsional oscilla¬ 
tion. This effect is akin to the seismic site effect lead¬ 
ing to the amplification of earthquake seismic waves 
in presence of certain geological conditions. In other 
words, considering a constant density one underesti¬ 
mates the matter displacement. 



Fig. 6.— Frequency of the f-modes with t = 1 as a 
function of a for different number of nodes: n — 1 
mode solid blue, n = 2 mode dashed red, n = 3 mode 
dotted green. These results have been obtained con¬ 
sidering the CCSC crust layer homogenous and the ra¬ 
dial dependence of the shear modulus and of the matter 
density of the ionic crust layer. 

As in the case of an homogenous ionic crust consid¬ 
ered in the previous section, the amplitude of the oscil¬ 
lations is very large and indicates that a crust cracking 
may occur. 

The shear strain corresponding to the ionic crust 
oscillation is reported in the upper panel of Fig. [8] as 
a function of the radius. Comparing Fig. [5] and the 
upper panel of Fig. [8] one can see that a new feature 
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Fig. 7.— Top: Amplitude of the oscillation of the fun¬ 
damental mode for a = 0.4. Both crust layers oscil¬ 
late with a comparable amplitude. Bottom: Ampli¬ 
tude of the oscillation of the fundamental mode for 
a = 0.8. The oscillation is basically segregated to the 
ionic crust. 


has appeared. The shear strain does not have a max¬ 
imum at the CCSC crust - ionic crust boundary. The 
reason is that the shear modulus of the ionic crust de¬ 
creases moving toward the surface and it is therefore 
easily displaced by a torsional oscillation. Flowever, 
the shear strain at the surface of the star must van¬ 
ish, because it corresponds to the no-traction boundary 
condition, thus a maximum is produced close to the 
star surface. 

This characteristic behavior of the shear strain is 
present also for the CCSC oscillations and the max¬ 
imum is located where the inhomogeneous term in 
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Fig. 8.— Shear strain produced by the \t\ oscillation 
when the inhomogeneous proprieties of the ionic crust 
are taken into account. Top: Shear strain across the 
ionic-crust layer for a = 0.8. Bottom: Maximum 
value of the shear strain produced for different values 
of the structure parameter a. 


Eq. © is maximum. In the bottom panel of Fig. [8] we 
show the maximum value obtained for the strain as a 
function of the parameter a. From this picture is clear 
that for any value of a the shear strain is large, larger 
than 10 -4 , possibly leading to crust cracking. We will 
discuss more in detail this issue in Sec© 

4.3. Temperature and magnetic field effects 

So far we have neglected the effect of the tempera¬ 
ture and of the compact star magnetic field. In this sec¬ 
tion we discuss the range of validity of the presented 
analysis; in particular we estimate the range of temper¬ 


ature and magnetic field for which the presented anal¬ 
ysis is approximately applicable. 


Regarding the temperature effects, it is known that 
the shear modulus decreases with increasing temper¬ 
ature, vanishing at the melting temperature. A typi¬ 
cal compact star surface temperature is of about 10 5 — 
10 6 K, but a larger te mperatu re is at tained in the inte¬ 
rior, see for example Shapiro & Teukolskv (1983). In 
any case, the typical energy scale of the CCSC crust is 
of order tens of Me V, thus the temperature effect in the 
CCSC crust layer can be neglected. On the other hand, 
a fraction of the ionic crust is believed to be strongly 
affected by the temperature. We will estimate the tem¬ 
perature at which the present analysis is valid consider¬ 
ing the amount of the ionic crust that is liquid, forming 
the compact star ocean. The existence of the ocean is 
due to the fact that the surface temperature of compact 
stars is larger than the melting temperature of some 
ionic crust chemical element. The transition from solid 
to liquid can be semiquantitatively determined by the 
ratio between the typical Coulomb energy and the ther¬ 
mal energy 

Z(r) 2 e 2 


r = 


r)T{r) 


(23) 


with the liquid/solid transition taking place at F = 
175. Therefore, the radial position of the liquid/solid 
transition point in the ionic crust depends on the local 
value of the temperature. Assuming a constant temper¬ 
ature in the ionic crust (which should be a very good 
thermal conductor), we can easily estimate the radial 
point corresponding to the solid/liquid phase transition 
in the ionic crust, see the upper panel of Fig. [9] For 
T = 3 x 10 7 K, corresponding to the dashed red line 
in the upper panel of Fig© only a small fraction of the 
crust is liquid. The value of the star density at which 
the transition takes place is p ~ 10'g/cm 3 . Therefore, 
in the present analysis we have assumed that the crust 
has at most a temperature of the order of 10' K. For 
larger temperatures, indeed, say for T = 3 x 10 8 K 
corresponding to the solid blue line in the upper panel 
of Fig© most of the ionic crust is liquid. 

Regarding the effect of the magnetic field, in the 
present model it is two-fold. There is a local effect 
on the net charge distribution located at r ~ R q and a 
standard magnetohydrodynamic (MHD) effect in bulk. 
The local effect corresponds to a deformation of the 
charge distribution at R q , and can be properly absorbed 
in the boundary conditions. Thus, it leads to a shift of 
the quantized frequencies. Unless one considers mag¬ 
netic fields B 10 14 G this shift is very small be- 
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Fig. 9.— Top: Coulomb parameter for two different 
values of the temperature. The horizontal line corre¬ 
sponds to T = 175, which is the value that separates 
the liquid and solid phases. Bottom: shear and Alfven 
velocities. The v a ,i curve is obtained for B = 10 13 
G, the va . 2 curve is obtained for B = 10 14 G and the 
curve u Qj 3 is obtained for B = 10 15 G. 


cause the electric field energy density in the electron 


density : 

nri986t) . 


layer is enormous (lAlcock et al.ll 1986T) . The bulk mag¬ 
netic field produces a Lorentz restoring force which 
competes with the shear stress. If the Lorentz restoring 
force is bigger than the shear stress the associated os¬ 
cillations are Alfven waves. In particular, if the Alfven 
velocity is much larger than the shear velocity, then 
MHD effects are dominant. In the bottom panel of 
Fig. [9] we compare the shear velocity with the Alfven 
velocity for three different values of the magnetic field. 
For B < 10 13 G the shear velocity dominates and the 


presented analysis applies. Note that Alfven waves are 
present even if the crust melts. 

5. Conclusions 

We have considered a model of nonbare strange star 
comprising color superconducting quark matter sur¬ 
mounted by a standard nuclear matter crust. The in¬ 
ternal part of the quarksphere is in the CFL phase; the 
external part of the quarksphere is in the CCSC phase 
and is separated from the ionic crust by a few hundred 
Fermi thick electron layer. We have determined the 
background configuration solving the pertinent TOV 
equation considering a simple parameterization of the 
EoS of quark matter and a realistic EoS for the descrip¬ 
tion of the ionic crust. We have considered one of the 
possible stellar structures in the sequence of configu¬ 
rations that solve the TOV equations. Similar results 
can be obtained considering strange stars with differ¬ 
ent masses and radii. 

Both the CCSC and the ionic crust are rigid and 
we have studied the torsional oscillations supported 
by these two electromagnetic ally coupled crusts. We 
have classified the torsional oscillations as CCSC crust 
oscillations and ionic crust oscillations, depending on 
which part of the structure has the largest displace¬ 
ment. CCSC crust oscillations are completely negli¬ 
gible if the CCSC crust is thin, say less than ~ 2 km, 
and only ionic crust oscillations are relevant. We have 
considered both a simplified model in which the two 
crusts are homogeneous and a more realistic model in 
which the ionic crust density and shear modulus are 
radial dependent. In any case we have assumed that 
the CCSC crust is homogeneous, because the density 
in the quarksphere changes by less than a factor 2. 
We have obtained that the oscillation frequencies of 
the £ = 1 modes are order of 10 kHz, and are not 
very sensitive to the extension of the CCSC crust. The 
£ = 1 modes correspond to oscillatory twists of the 
crust and do not conserve angular momentum, there¬ 
fore these modes are activated by events that transfer 
angular momentum to the strange star crust. A typ¬ 
ical event of this sort is a pulsar glitch associated to 
vortex drift from the CFL core to the CCSC crust. For 
that reason, we have assumed that the energy conveyed 
to the f-mode oscillation is of the order of the one of a 
Vela-like glitch. For definiteness we have assumed that 
all the energy of a Vela-like glitch triggers one single 
mode. The obtained deformation of the ionic crust is 
very large even considering CCSC crust oscillations. If 
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a fraction a of the Vela-like glitch energy is conveyed 
to the considered /-mode, then the amplitude and the 
shear deformation must be scaled by yja . If a ~ 1, 
the strain is such that it will possibly break the ionic 
crust layer. If the CCSC crust is sufficiently thin to 
segregate most of the oscillation in the ionic crust, a 
much smaller energy suffices to produce a large defor¬ 
mation of the ionic crust: even considering a ~ 10~ 3 
the shear strain on the ionic crust is of order 0.1. Crust 
cracking happens if the shear strain is larger than the 
breaking strain, s max , see Kittel ( 19761) for a discussion 


of the breaking strain in standard materials. The break¬ 
ing strain of the ionic crust is highly uncertain, indeed 
it is not known which is the microscopic mechanism 
responsible for the nonelastic response. The widely 
used range of values is between s max = 10~ 4 — 10 -2 , 
but values of 10“ 1 could be appropriate for perfect 
crystals without defects Kittell (1976!) as also shown by 
recent results obtained by molecular dynamics simula¬ 
tions of Coulomb crystals (Horowitz & Kadaul 2009). 
In order to properly asses where and when the crust 
breaks, it would be important to compute the breaking 
strain as a function of the radial distance. The result of 
our computation is, indeed, that the maximum strain 
produced by a it± oscillation is located at few tens of 
meters below the surface of the ionic crust. There¬ 
fore, this is the part of the ionic crust that will likely 
break during a glitch. We are not aware of any ob¬ 
servable that might be related to the breaking of the 
ionic crust at a specific radial distance . Possib ly, giant 
gamma-ray bursts, see I Thompson & Duncan (2001), 
and quasi periodi c os cillations, seellsraeLetak 2005ft 


Strohmav er & Wattsl d2005 l 2006 ): Glampedakis et al 
( 2006 ): Watts & Reddvl ( 20071) . might be put in rela¬ 
tion with such a phenomenon, but it would be inter¬ 
esting if some specific observable could help to dis¬ 
tinguish this breaking process from standard crustal 
breaks. 

Note that the maximum of the shear strain in 
the present model is much closer to the star sur¬ 
face than_jn_standard neutron stars, see for exam¬ 
ple McDermott et al. ( 1988b . Moreover, in standard 


neutron stars the energy of the t-modes is spread across 
the entire crust, extending for more than a km below 
the star surface. By contrast, in the considered model 
of nonbare strange stars, all the energy of the f-modes 
is conveyed in a layer few hundred meters thick, corre¬ 
sponding to the ionic crust at densities below neutron 
drip. Clearly, the effect of a reduced crust layer is to 
maximize the deformation. Considering that in our 


model we have taken the largest possible extension of 
the ionic crust layer, we have certainly underestimated 
the shear strain. 


As a final remark, note that if the f-oscillations are 
not triggered by angular momentum transfer to the 
strange star crust, then l= \ modes cannot be excited. 
Then, the low-lying oscillating mode is the %to mode. 
We find that the oscillation frequency of this mode 
is of the order of 5 kHz and it is weakl y dep ends on 
a. Similar results were reported in iLinl (120131) , where 
a nonbare strange star model with a CCSC core was 
considered. Assuming that a Vela-like glitch energy is 
conveyed to this mode, we find that the corresponding 
shear strain is suppressed with respect to the 1 = 1 
mode by about three orders of magnitude. We expect 
that similar results hold for higher values of i, meaning 
that assuming that an equal amount of kinetic energy 
is deposited to the modes with different value of i, the 
1=1 mode will produce the larger shear strain. 
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Fig. 10.— Torque acting of a system of two homogenous slabs in contact by a surface. We assume that the static 
friction between the two slabs is so large that the shear oscillations obey the no-slip boundary condition. 


A. Oscillations of two tightly bound rigid slabs 

Consider the system depicted in Fig. [lO] consisting of two slabs with a common contact surface. Suppose that a 
torque is applied to the structure. In standard conditions the two slabs would easily slide along the common surface 
with some kinetic friction. If the applied force is below a threshold value, the static friction prevents the relative 
motion and the whole system transversally oscillates around the equilibrium configuration. In common materials, the 
static friction is due to electrostatic Van der Waals forces and is extremely small. Suppose, however, that a strong 
electrostatic field is present along the contact surface, in such a way that the two surfaces are tightly bound. Let us 
consider the fictitious extreme case in which the binding force between the two surfaces is of the same order or large 
than the force between the atoms in the two slabs. Then, the force necessary to make the two surface slide is larger than 
the force needed to break the bounds in the bulk of the two slabs. The system of two slabs can be considered as one 
single nonhomogeneous body with a step-like matter density and shear modulus. In these conditions, the transverse 
displacement field, u, is a continuos function across the surface (no-slip boundary condition) 

ui(0-)=u 2 (0 + ), (Al) 

where we assumed that z = 0 corresponds to the contact surface. The no-traction boundary condition leads to 

vid z ui(Q~) = ^29 z m 2 (0 + ) , (A2) 


meaning that for unequal shear moduli, the derivative of the transverse displacement field is not continuous at the 
interface (a cuspid point). Note that the matter densities of the two slabs does not explicitly appear in this equation, 
indeed the shear moduli have dimension of a pressure. The matter density dependence is hidden in which in general 
is a function of the molecular binding force and density. The above no-traction boundary condition basically states 
that the transverse pressure on the contact surface must vanish. 

The frequencies of the quantized oscillations are obtained supplementing Eqs. (1A 1 l i and (IA2b with no-traction 
boundary conditions for the free surfaces, which for the considered planar geometry can be expressed as the Neumann 
boundary conditions d z ui(—Di) = 0 = d-wiiD-i)- In the linear approximation, a straightforward calculation gives 


(v\V2 + v 2 vi) sin 





(i/ 2 fi — V1V2) sin 





(A3) 
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where Vj = \]Vi/pi are the shear velocities. For v\v 2 v 2 Vi we obtain the approximate solution 


2 sin 




= 0 , 


(A4) 


showing that in this case there are two distinct frequency of oscillations, one determined by the mechanical property 
of slab 1 and one from the property of slab 2. Note that the frequency quantization 


^n 


nnvi 

“dT’ 


(A5) 


is the one that is obtained by imposing Dirichlet boundary conditions for a single slab with width D\ and shear velocity 
v\. On the other hand, the frequency quantization 


(2 n — 1)7TU2 
= 2 Lh 


(A6) 


is the one that is obtained by imposing Neumann boundary conditions considering a single slab with width D 2 and 
shear velocity v 2 . Note that we can recast the condition v-\ v 2 v 2 v-\ as 



(A7) 


and assuming that the slab 1 describes the CCSC crust and that the slab 2 describes the ionic crust, this inequality is 
certainly satisfied. Therefore, in the planar approximation, the oscillation can be divided in ionic crust oscillation with 
eigenfrequencies given in Eq. (lA5b . and in CCSC crust oscillation with eigenfrequencies given in Eq. ( |A6| i. As we 
have seen in Sec. 14. II this classification remains approximately valid in the spherically symmetric case. 
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